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Abstract 

In this second of the set of two papers on Lie symmetry analysis of a class of Lienard type 
equation of the form x + f{x)x + g{x) = 0, where over dot denotes differentiation with respect 
to time and f{x) and g{x) are smooth functions of their variables, we isolate the equations which 
possess maximal Lie point symmetries. It is well known that any second order nonlinear ordinary 
differential equation which admits eight parameter Lie point symmetries is linearizable to free 
particle equation through point transformation. As a consequence all the identified equations turn 
out to be linearizable. We also show that one can get maximal Lie point symmetries for the above 
Lienard equation only when f xx = (subscript denotes differentiation). In addition, we discuss the 
linearising transformations and solutions for all the nonlinear equations identified in this paper. 
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I. INTRODUCTION 



The present paper continues the work on the classification of Lie point symmetries de- 
scribed in the previous paper- of the second order ordinary differential equation of the 
Lienard type, 

A(x, x, x) = x + f(x)x + g(x) = 0, (1) 

where / and g are smooth functions of x and over-dot denotes differentiation with respect 
to t. In the first parti while solving the determining equations we have assumed the case in 
which one of the symmetry functions b(x) = and identified the corresponding equations. 
The question now naturally arises as to what are the invariant equations when b(x) ^ 0. In 
this part we present the answer to this question. 

On solving the determining equations with b(x) / we find that a class of equations 
possess eight parameter Lie-point symmetries. The forms of / and g which lead to the 
maximum number of symmetry generators are as follows (in increasing order of generality): 



(i) 


/ = 


0. 


9 


= Xix + A 2 , 


(ii) 


/ = 


k, 


9 


= X\x + A 2 , 


iii) 


/ = 


kx, 


9 


= -k 2 x 3 + Xix + A 2 , 
9 


iv) 


/ = 


kix + k 2 , 


9 


= \k 2 x 3 + \k x k 2 x 2 + Aix + A 2 
9 o 



where k, ki,k 2 , X\ and A 2 are constants. The first case includes the equation of motion of the 
free particle and simple harmonic oscillator. It is well known that both the equations admit 
eight symmetry generators' 1 ^ satisfying an sl(3,R) algebra. Similarly, the case (ii) which 
corresponds to the damped harmonic oscillator which has also been shown to admit sl(3, R) 
symmetry algebra 5 . From case (iii), by restricting the parameters Ai = A 2 = 0, one can get 
the modified Emden type equation which is the first nonlinear ODE shown by Mahomed 
and Leach to admit a rich set of Lie point symmetry generators^. Here, we show that the 
general modified Emden type equation with linear term and constant external forcing also 
admits eight point symmetry generators. We also report the class of equations which falls 
under case (iv) admitting eight symmetry generators and this may be considered as the 
most general form of ([I]). It is important to note that all these equations are linearizable 
under point transformations. The explicit forms of the infinitesimal symmetries and their 
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associated generators for the general cases are also reported for the first time, as far as the 
authors' knowledge goes. 

In all the above cases, it turns out that f xx = 0. By considering the next higher degree 
polynomial in the variable x for /, say quadratic, / = kix 2 + k 2 x + k%, so that f xx 7^ 0, we 
get 

g = ^klx 5 + j^k^x 4 + Aix 3 + X 2 x 2 + X 3 x + A 4 , 

which leads to one symmetry generator only. In general, we prove that for f xx ^ one 
obtains lesser parameter symmetry group only. 

The plan of the paper is as follows. In Sec. II we recall the derivation of the determining 
equations for the infinitesimal symmetries and describe a procedure to solve them. In Sec. Ill 
by demanding all the four symmetry functions to be non-zero, we construct systematically 
the class of functions / and g which admit eight symmetry generators. We start our analysis 
with / = 0, g = in Eq. (CQ), which lead to the equation of motion of the free particle, and 
derive all earlier known equations from (TjQ) which admit eight point symmetries. Proceeding 
further by assuming the function / to be a constant we deduce damped harmonic oscillator 
equation and its variant equations. We also present the infinitesimal symmetries and their 
associated generators for these equations. Next, in Sec. IV, by considering the function / to 
be linear in x, we deduce two new family of nonlinear ODEs which admit eight symmetry 
generators. In Sec. V, we prove from the symmetry determining equations that when f xx ^ 
only lesser parameter symmetries can exist. We also show the equivalence of our results with 
the Lie's linearization criterion. Finally in Sec. VI we provide a summary of our results. 
The solutions of the determining equations encountered in the symmetry analysis of new 
equations are presented in Appendices A and B. 

II. SYMMETRY DETERMINING EQUATION OF flU WITH b^O 

To make this paper self contained we recall the essential determining equations very 
briefly here. The invariance of Eq. (CQ) under one parameter Lie group of infinitesimal 
transformations, t — t + e£(i , x) + 0(e 2 ), x = x + er](t, x) + 0(e 2 ), e C 1, leads us to four 
determining equations for the functions £ and r\ (vide Eqs. (8)-(ll) in Paper I). The first two 
determining equations can be integrated straightforwardly to yield £ = a(t) + b(t)x, rj = 
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bx 2 — 263(x) + c(t)x + d(t) , where 3^ = F(x) = f Q f(x')dx' and 3 XX = f(x). The remaining 
two determining equations, after substituting the forms of £ and T], read 

{bx 2 - 263 + cx + d)f x + (d + 6x)/ + 36# - 46F + 36a; + 2c - a = 0, (3) 

and 

(bx 2 - 263 + cx + d)g x - (c - 2d - 26F)# - 263 

+(&z 2 - 263 + cx + d)f + 'bx 2 + cx + d = 0. (4) 

Since g is linear in (j3j) we can rewrite ([3]) to obtain 

g = \\-(bx 2 - 263 + cx + d)f x - (d + bx)f + 46F - 36x - 26 + a], b + 0. (5) 
36 

While deriving Eq. (jSj) we assumed that 6 7^ 0. On the other hand one may also consider the 
case 6 = 0. In fact this was the case of the Paper I. Now we solve the determining equations 
(JSJ) and (jlj) with 6 7^ and identify the functions that are invariant under eight parameter 
Lie point symmetries. 

III. MAXIMAL LIE POINT SYMMETRIES OF LIENARD TYPE SYSTEMS: 
LINEAR ODES 

To begin with we consider polynomial forms for / in the variable x with increasing 
degrees and derive the corresponding explicit forms of g through Eq. (jSJ). Substituting 
the expressions of / and g into Eq. (jlj) and solving the resultant equation we obtain the 
infinitesimal symmetries £ and rj which in turn leads to the classification of the class of 
second order ODEs which admit maximum number of symmetry generators. We find that 
the maximal number (eight) of symmetries exists only when f xx = 0. When f xx 7^ 0, lesser 
parameter Lie point symmetries only exist. 

A. Linear Undamped Systems: / = 

To start with we consider the simple form / = 0. Correspondingly from Eq. (jHj) we get 

g = \ 1 x + \ 2 , (6) 
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where Aj and A2 are arbitrary constants. Substituting Eq. (jSj) into (jlj) and solving it con- 
sistently we can find the infinitesimal symmetries. However, to classify the results system- 
atically we consider the following subcases, viz., (i) Ai = A 2 = 0, (ii) A 2 7^ 0, Ai = 0, 
(in) Ai ^ 0, A 2 = and (iv) A x ^ A 2 ^ 0. 



1. Free particle motion (Ai = A 2 = 0) 

When both Ai and A 2 are equal to zero, the function g also becomes zero and so Eq. ([1]) 
takes the form 

x = (7) 

which is obviously the equation of motion of a free particle. Substituting / = g — into 
Eq. (jlj) we obtain the equation b x 2 + cx + d = 0. Equating the coefficients of x\ i = 0, 1, 2, 
to zero separately one obtains 6 = 0, c = and d = 0. Solving the later two equations one 
obtains c = c\ + c 2 t and c? = d\ + <i 2 £, where d { s and d^s, z = 1, 2, are arbitrary constants. 
Similarly from Eq. ([3]), we obtain 6 = and a = 2c. Solving these two equations one obtains 
b = bi + b 2 t and a = a\ + a 2 t + c 2 t 2 , where b^s and a^s, z = 1,2, are arbitrary constants. 
Inserting the forms a(t), b(t), c(t) and d{t) in £ = a(t) + 6(t)x and 77 = 6a; 2 + c(t)x + one 
gets 

£ = a x + a 2 t + c 2 t 2 + (61 + b 2 t)x, t] = d x + d 2 t + (ci + c 2 t)x + b 2 x 2 . (8) 
The associated eight infinitesimal generators take the following form 

Xi= ^ X2=t ^ Xs=x ^ X4= H^ +X ^ 

9 / <9 <9 A 9 9 

X 5 = x— , X 6 = t t— + x— , ^7 = ^, ^8 = ^, (9) 
ax \ at ax / ax ax 



he symmetries and their generators (Q coincide exactly 



which satisfy an si (3, R) algebra, 
with the known ones given in Refs 



2. Free falling particle (A 2 , ^ 0, Ai = 0) 

If we choose Ai = and A 2 as arbitrary in Eq. ([6]), then from Eq. (CQ) we get 

x + A 2 = 0, (10) 



which is the equation of motion of a free falling particle. Now solving the determining Eq. 
(j3J) with / = and g = A 2 we obtain^ 

3 A 
f = ai + a 2 t + (c 2 + -\ 2 h)t 2 + y 6 2 t 3 + (61 + 

A 1 A 2 

77 = + d 2 t + y (ci - 2a 2 )t 2 - A 2 (A 2 6! + -c 2 )t 3 - -j-b 2 t 4 + ( Cl + c 2 t)x + 6 2 x 2 , (11) 

where a^fo^Cj and c?j, i — 1,2, are arbitrary constants. Here too, we have the following 
eight generators 2 - 

* = - Y -4- < i2 > 

which also form an sl(3,R) algebra. 



3. Free linear harmonic oscillator (Ai 7^ 0, A 2 = 0) 

We consider Ai as arbitrary and A 2 = so that Eq. (JTj) becomes 

i + Aia; = 0, (13) 

which is nothing but the equation of a linear harmonic oscillator (for Ai > 0). Solving 
Eq. (jl]) with / = and g = \\X we get the following well known infinitesimal symmetry 
transformations 2 

£ = ai + a 2 sin2at + 03 cos2at + {b\ sin at + 6 2 cos at)x, (14) 
Tj — d\ sin at + c? 2 cos at + (ci + a 2 a cos 2at — a^a sin 2at)x + a(b± cos at — 6 2 sin at)x 2 , 

for Eq. ( fl~3~I) where a 2 = Ai and a,, 6j, cfj, z = 1, 2, 3, j = 1,2 and ci are arbitrary constants. 
The associated generators forming the sl(3,R) algebra are 2 

d d d d d 

Xi = — , X 2 = sin2at— + axcos2at— , X 3 = cos2at— axsin2at— , 

at Ot ox ot ox 

d d d d 

X 4 = x(sin at— + ax cos at) — ,X 5 = x(cos at- axsinat) — , 

ot ox ot ox 

Odd 
X e = x—, X 7 = sin at— , X 8 = cos at— • (15) 
ax ox ox 
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On the other hand in the case Ai < 0, Eq. (|T3|) becomes the repulsive harmonic oscillator 
which is invariant under the following infinitesimal symmetries 

f = ai + a 2 e 2at + a 3 e~ 2at + (b ie at + b 2 e~ at )x, 

v = d ie at + d 2 e~ at + (ci + a{a 2 e 2at - a 3 e- 2at ))x + a(b x e at - b 2 e" at )x 2 , (16) 
The associated sl(3,R) symmetry generators are 

A" 5 = le-"'^- - ax^-), A' e = e°'|-, A' 7 = e -'|-, X, = x%-. (17) 
at ox ox ox ox 

4- Displaced linear harmonic oscillator (Xi ^ X 2 / 0) 

Finally, we consider both Ax and A2 to be arbitrary (not equal to zero) and we have an 
equation of the form 

x + Aix + A 2 = 0, (18) 



which corresponds to a displaced simple harmonic motion. Eq. (jTBj) can be transformed 
to (1131) by a simple transformation, x — > x'= x + (^) and consequently the infinitesimal 
symmetries of the former can be derived from that of Eq. (jl4jl . 



B. Linear Damped Systems: / = constant = k 

In the previous subsection, III A, we considered the case / = and derived the form of 
g. Now we fix / = constant = k and deduce the associated form of g. Substituting / = k in 
Eq. ([5]), we again obtain g as a linear function in x, that is 

= A x :r + A 2 , (19) 

where Ai and X 2 are arbitrary constants. Consequently one can solve Eq. (j3J) with (|T9l) to 
obtain the infinitesimal symmetries. However, to relate the present results with those of the 
literature, we consider the following four cases separately: (z) Ai = A 2 = 0, (ii) Ai 7^ 0, A 2 = 
0, (Hi) A 2 ^ 0, Ai = 0, and (iv) Ai^ A 2 ^ 0. 
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1. Free particle in a viscous medium (Ai = A2 = 0) 

Let Ai = A2 = 0. Then from Eq. (fT9|) we obtain g = 0. Substituting f = k and g = in 
Eq. (jlj, we get 

x + kx = 0, (20) 

which corresponds to the equation of a free particle in a viscous medium. Solving Eq. (J3J) 
with f = k and g = we get 

i = ^ + T eU + ¥^ ki + ^ + b i eU)x - 

!) = ((,+ ^e-* 1 + (c, - - H^ 2 , (21) 

where a i; fej, q and c?j, z = 1,2, are arbitrary constants. The corresponding generators 
satisfying the sl(3, R) algebra are 8 

X - 9 X - l p ktd X-r( d kr 9 \ X - X P ktd 

x * = 4' x * = V kt ( \% - 4) k kt f ( 22 ) 

ox k \k at ox J ox k ox 

2. Damped linear harmonic oscillator (Ai 7^ 0, A2 = 0) 

Considering Ai as arbitrary and A2 = in Eq. (Tl9l) . we obtain from (0Q), the equation of 
a damped harmonic oscillator, 

x + kx + X 1 x = 0. (23) 

The invariance property of this equation has been discussed in detail in Reffl Following 
our procedure we find Eq. (1231) to be invariant under the infinitesimal transformation, 

£ = ai + a 2 sin 2(3t — a 3 cos 2(3t + e at {pi cos j3t + b 2 sin/?i)x, 
?7 = [&3 + a 2(/3 cos 2/3t — a sin 2/3t) + a^{(3 sm2(3t — a cos 2/%)] x 
- [6ie at (a cos (3t + (3 sin (3t) + b 2 e at (a sin /ft - (3 cos /ft)] x 2 
+e" a * (di cos /ft + d 2 sin /ft ) , (24) 

where a = |, /3 = |(4Ai — A; 2 )5 and a i; foj, dj, i = 1, 2, 3, j = 1, 2, are arbitrary constants. 



S 



The sl{3, R) symmetry generators^ are 

d d d 

X x = — , X 2 = sin2/3t— + x(/3cos2/% - a sin 2/%) — . 
dt dt ox 

d d 
X 3 = -cos 2/%— + x(3sm28t - acos23t) — , 
at ox 

( d d 

X 4 = xe at ( cos Bt— x(a cos 6t + 3 sin /3i) — 

\ at ox 

( d d 
X 5 = xe at I sin 3t— + x(3 cos 3t — a sin /3t) — 
\ dt Ox 

r\ r\ r\ 

X 6 = x—, X 7 = e- at cos/?t— , X 8 = e-^sin/^t— . (25) 
ax ax ax 

In the above derivation, 4Ai < k 2 . Then for the choice 4Ai = A; 2 , we have the following 
infinitesimal symmetries: 

£= ai + aat + aa^ + ^ + ^a + Me - ^)^ (26) 
^ = - 0(&2 + M) - M zV** + ^ + (a 3 - |a 2 )t - | a 3^ x + (d l + d 2 t)e^ 1 . 

The respective infinitesimal vector fields forming an sl(3, R) algebra are 

d d k d J k . d 

Xi = — , X 2 = t- tx— , X 3 = t 2 — + xt(l 1) — , 

dV dt 2 dx dt v 2 1 dx 

ktd v _* t (d 3, d\ „ _» t A9 3, 9 

Finally the choice 4Ai > k 2 leads to the following infinitesimal symmetries: 

£ = a x + ^ - ^e"^ + fae^' + 6 2 e( Q -«*)x, 

rj = ({p-a)b 1 e^ t -{3 + a)b 2 e^ t )x 2 + d 1 e^ a+ ^ t 

+ ( Cl + 1(1 - j)a 2 e^ + 1(1 + ^W"*) x + d*-^, (28) 

where a = |, /3 = |(/c 2 — 4Ai)s. The associated vector fields forming an sZ(3, i?) algebra are 

Xl - at' X2 " e U^ + 2 (1_ y )x ^ 



= *,-(«+£)* JL 

O ^ Q 

ax 
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It is well known that equation (1231) admits a time dependent Hamiltonian of the form 
H = §^e~ kt + ^x 2 e kt . However, very recently we have proved that the system (123]) also 
admits a conservative Hamiltonian description^ for all the values of k and Ai. The explicit 
form of the Hamiltonian can be deduced from equation (71) of paper I by fixing a = k, 
(3 — Ai and q — 1 in it. 

3. Falling particle in a viscous medium (A2 7^ 0, Ai = 0) 

Choosing Ai = and A2 to be arbitrary in Eq. ( 1T91) . we obtain from Eq. (JT]) 

x + kx + A 2 = 0, (30) 

which is the equation of a falling particle in a viscous medium. Eq. ( 1301) can be transformed 
into (120]) through the transformation x — * x' = x+ (4r)i, and so its infinitesimal symmetries 
and generators can be derived from fpZTT) and (1221) by appropriately replacing the variable x. 

4- Displaced damped harmonic oscillator (Ai 7^ A2 / 0) 

Let Ai ^ A 2 ^ in Eq. (US}. Then Eq. © becomes 

x + A;x + Aix + A 2 = 0. (31) 

Here too one can transform Eq. (|3T|) to the damped harmonic oscillator (1231) through the 
transformation x — > x' = x + ^ and so the invariance properties can be analyzed from the 
infinitesimal symmetries of the damped harmonic oscillator. 

IV. MAXIMAL LIE POINT SYMMETRIES OF LIENARD TYPE SYSTEMS: 
NONLINEAR ODES 

A. Modified Emden Equations: f = kx 

From the above study, we find that the choices / = and / = constant lead to linear 
ODEs only. Thus to generate nonlinear ODEs which admit maximum number of symmetry 
generators one needs to consider linear/higher degree polynomial/nonpolynomial forms for 
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the function f(x). To begin with we consider / to be a linear function of x, say / = kx, 
where k is a constant. As a consequence we get (from Eq. (jSJ)) 

g = IfcV + \ lX + A 2 , (32) 
9 

where Ai and X 2 are arbitrary constants. Substituting the explicit forms of / and g into 
Eq. (j4j) and solving it consistently we obtain the symmetry generators of Eq. (QQ). For 
this purpose, we consider the following four cases, (z) Ai = A2 = 0, (zz) Ai 7^ 0, A2 = 0, 
(Hi) A2 7^ 0, Ai = and (iv) Ai 7^ A2 7^ 0, separately. 



i. Modified Emden equation: Ai = A2 = 



For Ai = A2 = 0, we obtain from Eq. (CQ) the following nonlinear ODE 

k 2 

x + kxx + ~^~ a;3 = 0? 



(33) 



whose linearization, invariance and integrability properties have been widely discussed in 
the recent literature^ 1 ^ 1 ^^ 3 - 1 ^. From our above analysis, Eq. (133]) is found to be invariant 
under the Lie point symmetries 



k k 
£ = a x + a 2 t + (c 2 + -d^t 2 - -d 2 t 3 + 

2 o 



61 + b 2 t - -(ci + a 2 )t 2 
o 



k k k 2 

— (C2 + -dx)t 3 H d 2 t 4 

6 V 3 ; 36 



k 



x, 



77 = c?i + d 2 t + (ci + c 2 t - -d 2 t 2 )x + 



k 



k 



k k 2 
b 2 ~ |(ci + a 2 )t - '-(c 2 + -d,)t 2 + —d 2 t 3 



x 



k 
3 



h + b 2 t - |(ci + a 2 )t 2 - ^(c 2 + ^x)t 3 + ^-rf 2 t 4 
o o 3 3d 



x 



(34) 



where a^, 6j, q and g?j, z = 1, 2, are arbitrary constants. These infinitesimal symmetries lead 
to the symmetry generators 

9 /„ k . d k . k . 9 9 k n d 

X x = —, X 2 = t(l art)- x 2 t(l xt) — , X 3 = x- a; 3 — , 

dt' V 6 ! dt 3 V 6 ; &r' 9t 3 



X 4 
X 6 



9 9/ „ N d 
xt—+x 2 (l--xt)—, X 5 
at 3 ox 



--xt 2 — + x(l xt H x 2 r) — . 

6 v 3 18 ; <9x' 

O / - k \ d /. k k ry n Q 
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x 7 = ^ 2 (i - ±xt)?- + (1 - -tV + -tV)^-, 

7 2 v 9 ; <9t v 6 54 ; &r 

X 8 = --t 3 (l - -zt)-J + til - -xt + —x 2 t 2 - —x 3 t 3 )-^-, (35) 

8 6 v 6 J dt K 2 9 108 ; <9x v ; 

which also satisfy an sl(3, R) algebra 6 . Eq. (133]) has the general solution 

= ^+y t ' + / 2 - < 36) 

where Ji and I2 are two integrals of motion with the explicit forms 

h = ~t + 1 -^—, J + (37) 

2. Modified Emden equation with linear term (Ai ^ 0, A2 = 0) 

Choosing Ai as arbitrary and A 2 = in Eq. (1321) . we obtain 

2 = —x 3 + Xtx. (38) 
y 

As a result, we have a nonlinear ODE of the form 

k 2 

x + kxx + —x 3 + Aix = 0, Ai ^ 0. (39) 
9 

Using the forms of / and g in Eq. (T4|), we obtain the following set of coupled ODEs for the 
four arbitrary functions a, b, c and d: 

36 + 3Ai6 + ka + kc = 0, 2c-d + £;d = 0, c + 2Aid + kd = 0, 

3 b + 3Ai6 + 3A;c + k 2 d = 0, d + Aid = 0. (40) 

Note that the fourth equation in (14131) can be obtained from first and second equations and 
one has effectively only four equations for the four unknowns. Solving equations ( 1401) we 
find that eight parameter Lie point symmetry groups exist for both the choices (i) X± > 
and (ii) Ai < 0. In the first case the infinitesimal symmetries read 



£ = di sin 2at ^- cos 2at -d\ sin at —^di cos at + (61 sin at 

2a 2 2a 2 3a 2 3a 2 V 



k k kc^ 

+6 2 cos at H -CI3 cos 2at H -a 2 sin 2at -)x, 

6a 2 6a 2 3a 2 

(ci H cos 2at H sin 2at H aj cos at ofo sin at)x 

2 2 3a 3a 

+ (61 a cos at — 62a sin at 03 sin 2at H a 2 cos 2at)x 2 

3a 3a 



1 k k kc~^ 

■-(61 sin at + b 2 cos at H -a 3 cos 2at H -a 2 sin 2at -)x 3 

3 6a 2 6a 2 3a 2 

+di sin at + d 2 cos at, (41) 
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where a%, bi, Cj, di, i = 1, 2, j = 1, 2, 3, are arbitrary constants and a = yAi- 
The corresponding symmetry generators are 



X 1 

x 3 

x 4 
x 5 
x, 
x« 



9 „ , ■ o k n , d . n 2k . n 
— , A 2 = sin2at + — xcos2at)— - + x(a cos 2at — — zsin2at 
at 3a at 3 

k 2 2 n v <9 

a; cos 2at) — , 

9a ax 

fc .9 . 2k 

(cos 2at x sin 2cm x(a sin 2at H x cos 2at 

v 3a 'dt y 3 

/c 2 2 . d 
— —x sm2at)— — , 
9a ax 

d k d 

x sin at— + x 2 ( a cos at x sin at) — — , 

at v 3 ' dx 

° it k .8 

x cos at— x (asm at H — x cos at) — , 

at 3 ax 

9 ,3a 2 x 1 o. 9 fc 8 , . k .a 

x— - (— h -x )— , A 7 = -— — smat- + (sin at + — xcosat) — , 

at k 3 ax 3a 2 at 3a ax 

k d k d 
cos at— + (cos at x sin at) — 



(42) 



3a 2 dt 3a ax 

The vector fields (I42p can be shown to form an sl(3,R) algebra. For the case (ii) the 
infinitesimal symmetries involve exponential functions of t, that is 



k 



k 



k 



£ = ai + a 2 e 2at + a 3 e~ 2at + —;d x e at + —:d 2 e~ at + [— c x + he at + b 2 e 



-at 



3a 



3a 2 



3a 
1 



[a 2 e 2at - a 3 e- 2at 



))*, 



*p- 2 + he at + b 2 e~ at - A(a 2 c 2Qi - a 3 e~ 2at ) 
I3a z 3a 



L 3a 2 



x + [afeie - ab 2 e 



(43) 



-ai 



2k 



(a 2 e 



2at 



+a 3 e 2at )]x 2 + [ci + a(a 2 e 



**t _ a3e -2at. 



— (dte * - d 2 e- at )]x + die * + d 2 e . 
3a 



-at 



where a^, C\,bj, dj, i — 1, 2, 3, j = 1, 2, are arbitrary constants. 

The corresponding symmetry generators, which also satisfy an sl(3,R) algebra, are 



X 1 

x 3 



x 6 

X* 



8_ 
d~f 



X 2 = e 



2o/ 



1 - 



k \ 8 



3a X J dt 
k 



-2at 



k \ 8 



+ 



k g 2/c 2 
9a 3 



8_ 

dx 



3a J dt \9a 



2k 



1 + —x I —x + —x + ax 7— 



a 



X * = X di 



d tk 



3a 2 \ d 



r 



-X 



-at 



-at 



X 



d_ 

dt 



k J dx 
k 

3 i 2 

3 



d_ 

dx 



d 



< 3a 2 \ d i 
dt V k J dx 



X R 



dx. 

at Q . 



dt 



—x — ax 



X 7 = e 



at 



8 _ / _ 3a 2 \ d 
.dt V k J dx 



2 \dj 
) dx 



(44) 
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The nonlinear oscillator equation (139]) admits several unusual dynamical properties. In 
fact the dissipative system fl39l) admits a time independent first integral and a conservative 
Hamiltonian description 15 for all values of k and A, that is 



where 



_dL _27AV 1 A 3A 



For A > 0, the system (1391) admits explicit amplitude independent sinusoidal periodic 
solution of the form 

x(t) = r — i — ; — '- — -, < A< -, UJ = V\ 47 

and when A < 0, the solution turn out to be a front like one^ 5 -, that is 

3/3(7,^ - 1) 

W t/i/ 2 e* + t(l + /ief)' 11 
where /3 = ^/\X\ and A, 5, Ii, I 2 are constants. 



5. Modified Emden equation with constant external forcing (A2 7^ 0, Ai = 0) 

In Eq. (13*21) let us put Ai = and A 2 = arbitrary so that Eq. (""J becomes 

k 2 

x + kxx + —x 3 + A 2 = 0. (49) 
9 

Substituting f = kx and g = yx 3 + A 2 in (TJj), we obtain the following determining 
equations for the arbitrary functions a, 6, c and d: 

36 + ka + kc = 0, 2c — a + 3A 2 6 + fed = 0, c + fed = 0, 

d- A 2 c + 2A 2 a = 0. (50) 

Solving them consistently one obtains 

a(t) = ai - A;(mi6 2 e Pl ' + m 2 6 3 e P2 * + m 2 6 4 e P3 ' + mi6 5 e P4 ' + mi6 6 e Ps * + m 2 6 7 e P6 '), 
6(0 = h + ^ + ^e P2t + ^e P3 < + ^e P4 < + + *Ie»«, 

Pi P2 P3 Pi P5 P6 

c{t) = ~k(m 2 p 1 b 2 e Plt + m 1 p 2 b 3 e P2t + m 1 p 3 b i e p ' it + m 2 p^b 5 e Pit + m 2 p 5 6 6 e P5 ' + m 1 p 6 b 7 e mt ), 
3A 

d(t) = r^h + n6 2 e P1 ' + r 2 6 3 e P2 * + r 3 6 4 e P3 * + r 4 6 5 e P4 ' + r 5 6 6 e P5 * + r 6 6 7 e P6 *. (51) 
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The above derivation is presented in Appendix A. Here ai and bi, i = 1, . . . , 7, are integration 
constants which are also the eight symmetry parameters. The constants Pi and 7,, where 
i = 1, . . . , 6, are introduced for simplicity which depend on the system parameters and can 
be fixed from the relations (IA15[) . and (IA20I) . Finally, the constants m\ and m 2 are defined 
in the expression (1A19I) . 

Eq. (|5~T!) provides us the infinitesimal symmetries of the form 



£ = ai - k(mib 2 e Plt + m 2 b 3 e P2t + m 2 6 4 e P3t + mib 5 e P4t + mib & e 



P3t 



m 2 b-je 



+ [b 1+ ***** + + b -e 

Pi P2 P3 



Pi P5 P6 



x. 



(6 2 e pi * + b 3 e P2t + b 4 e P3t + b 5 e P4t + b 6 e P5t + b 7 e 
+ h e p 2 t + h e p>t + ^ + h e pst + ^ 



jfc 



77 ( 61 + -e 
Pi 



.r 



{m 2 p\b 2 e 



Pit 



P2 P3 Pi P5 P6 

+mip 2 6 3 e P2< + m 1 p 3 6 4 e P3 * + m 2 p i b^e P4t + m 2 p 5 b e e P5t + m 1 p e bje P6t ) x 
3A 2 



-61 + nb 2 e Plt + r 2 b 3 e P2t + r 3 6 4 e P3 * + r 4 b 5 e P4t + r 5 b 6 e Pst + r 6 b 7 e 



Pet 



The associated infinitesimal generators are 



X 1 
X 3 

x 4 
x 5 
x 6 
x 7 

X* 



d_ 



d 



— , X 2 = x— - -A 2 + -x 6 — 



oPlt 



oPrf 



oP3t 



oPit 



oP5t 



oP6t 



dt 



—krrii + 



—km 2 + 



—km 2 + 



-km\ + 



—km\ + 



—km 2 + 



k 







d_ 
Pi J dt 



dx 



+ I T\ — m 2 p\X + x — 



X 



d 



p 2 J dt 
x \ d 



+ r 2 — m\p 2 x + x — 



p 3 J dt 



+ r 3 — m\p 3 x + x — 



x 



d 



Pa J dt 
x \ d 



+ r 4 — m 2 p4X + x — 



PbJ dt 
x \ d 



+ r 5 — m 2 p 5 x + x - 



p 6 J dt 



+ I r 6 — mipQX + x — 



3pi 

A 

3p 2 
fc 

jfc 

3p 4 

A 

3p 5 

jfc 

3p 6 



dx 

d_ 

dx 

d_ 

dx 

d_ 

dx 

d_ 

dx 

d_ 

dx 



The general solution of Eq. (|49j) can be deduced explicitly following Ref. 

18a 2 (1 - J ie ±iQv/ 3*) 
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as 



x(t) 



—3a 1 
~k~ + 3k 



3a{l - I ie ±ia ^) T 2iay/3IJ 2 e 



± iaV3(l + 1^™^) 



(52) 



(53) 



(54) 
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where a is solution of the equation a 3 + X\a — |A 2 = and I\ and J 2 are two integrals of 
motion whose explicit forms have also been reported in Refil6. 



4- Modified Emden equation with linear term and constant external forcing (Ai 7= 

0, A 2 ^ 0) 

Finally, we choose Ai 7^ A 2 7^ in Eq. fl32|) . so that 

k 2 

g = —x 3 + Xxx + A 2 . 
y 

For this case Eq. (CQ) becomes 

x + kxx + -k 2 x 3 + Aia; + A 2 = 0. (55) 
9 

The determining equations for the infinitesimal symmetries become 

36 + kd + 3Ai6 + kc = 0, 2c - a + 3X 2 b + kd = 0, 

c + 2Aid + kd = 0, , c/-A 2 c + 2A 2 a + Aid = 0. (56) 

Even though Eq. (1561) looks similar to (1501) . the appearance of additional terms makes it 
tedious to solve. We solve Eq. (J5B]1 with the assumption that fc, Ai, A 2 7^ to obtain 
2A 2 

a W = fl i - 7TT 6 i + Pihe ait + p2b 3 e~ ait + b 4 e a2t {(3 3 cos a 3 t + (3 A sin a 3 t) 
k 2 X 2 

+b 5 e a2t (/3 3 sin a 3 t — /3 4 cos a 3 £) + b 6 e~ a2t (f3 5 cos a 3 t + /3 4 sin a 3 t) 

+i>7e _a2 * (/?4 cos a 3 t — /3 5 sin a 3 t) , 

6(t) = 61 H e ait e~ ait + — ^-(a 2 cosa 3 t + a 3 sino; 3 t) 

ai ai (a| + a 3 ) 

+ 7-^5 or (a 2 sin a 3 t — a 3 cos a 3 t) + -p—= ~- (a 3 sin a 3 t — a 2 cos a 3 i) 

(a| + a|) (« 2 + « 3 ) 

+ 7—9 97 ( a 2 sin a 3 t + a 3 cos a 3 t) , 

(«2 + "3) 
3A 

c(t) = -ifei + 7i&2e aii + l2b 3 e~ ait + 6 4 e a2 *(7s cos a 3 t + 74 sin a 3 t) 

+6 5 e a2 *(7 3 sin a 3 t — 74 cos a 3 t) + b 6 e~ a2t (75 cos a 3 t + 7 6 sin a 3 i) 

+6 7 e~ Q2 *(7 6 cosa 3 t - 75 sin« 3 t), 
3A 

d(t) = r^&i + 5i6 2 e Ql * + 5 2 b 3 e~ ait + b A e a2 H8 3 cos a 3 t + 5 4 sin a 3 t) 

k 

+b 5 e° 2t (6 3 sin a 3 i — 5 4 cosa: 3 t) + 6 6 e _Q!2 '(55 cosa 3 t + 5 6 sina 3 t) 

+b 7 e~ a2t (S 6 cosa 3 t — 5 5 sin a 3 t). (57) 
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Here a\ and bi, i = 1,...,7, are the arbitrary constants arising from the integration of 
Eq. (j56p which fix the eight symmetry generators. The other constants are included for 
simplicity sake. A detailed analysis on the above derivation and the relation between these 
constants and system parameters are given in Appendix B. 
The infinitesimal symmetries turn out to be 

2\j 



6 = 01 



k 2 \- 



h + Pib 2 e ait + (3 2 b 3 e~ ait + b A e a2 \f3 3 cosa 3 t + (3 4 sma 3 t) + b 5 e a2t (/3 3 sma 3 t 



— /3 4 cos a 3 t) + b 6 e ° 2 ((3 5 cos a 3 t + /3 4 sin a 3 t) + b 7 e ° 2 (/3 4 cos a 3 t — (3 5 sin a 3 t) 



-\-b-je 012 (/?4 cos a 3 t — (3 5 sin a 3 t) + 



0L\ Gil («2 + a 3 ) 



+a 3 sina 3 t) + ~- [a 2 sin a 3 t — a 3 cos a 3 t) + — -= „— [a 3 sin a 3 t — a 2 cos a 3 t) 



(«2 + «§) 

(a 2 sin a 3 t + a 3 cos a 3 t) 



(«2 + «i) 



X. 



3 



62 
a 1 



b A e a2t b 5 e a2t 

^- (ct 2 cos a 3 t + a 3 sin a 3 t) + 7—= ^- (a 2 sin a 3 t 

(«2 + a|) («i + ag) 



-a 3 cos a 3 t) + 

(«2 + «3 



,, v a 3 sin a 3 t — a 2 cos a 3 t) + 2 — ^- (a 2 sin a 3 t 



+a 3 cosa 3 t)] x 3 + 



x 2 + 



«2 + «3) 

b 2 e ait + b 3 e' ait + e a2t (b 4 cos a 3 t + b 5 sin a 3 t) + e~° 2 * (6 6 cos a 3 t 
3Ai 



+ 7i&2e ai ' + 7 2 6 3 e Ql * + & 4 e a2< (7 3 cos a 3 t + 74 sin a 3 t) 



—b 7 sin a 3 t) 

+6 5 e Q2 *(7 3 sin a 3 t — 74 cos a 3 t) + &6 e ~ Q2 *(75 cosa 3 t + 7 6 sin a 3 t) + b 7 e~ a2t {^Q cos a 3 t 
-75 sin a 3 t) 

+br,e a2t (5 3 sin a 3 t — £4 cos a 3 t) + &6e~ Q2 *(<55 cos a 3 t + 5 6 sin a 3 t) 

+bje^ a2t (S 6 cos a 3 t — 5 5 sin a 3 t). (58) 



x - -r^fci + 5i6 2 e aii + 5 2 b 3 e~ ait + b 4 e a2t (5 3 cos a 3 t + <5 4 sin a 3 i) 



The associated generators constituting a sl(3,R) algebra are 
d ( 2 At \ <J / 3A 2 3Ai 



X\ — — , A 2 



2A? \ 9 

X — T77TTT I TT7 + 



k 2 \\) dt 



k 3 \ d 
k k 3 / dx : 



X 3 

x 4 



x \ d 



k 



& + — ) 7T+ + Ui+7i^ + ^ 2 + q ... 



9 



x \ d 



k 



@2 ~ — ) + [82 + 12X + X 2 - X' , 

ol\) at V dai / ax 







-a-\ t 
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x R 



f3 3 cos a 3 t + /3 4 sin a 3 t 



1 . , U 

p-s- — rr- (a 2 cos a 3 t + a 3 sin a 3 t)x — 



+ (#3 cos a^t + 5 4 sin a 3 t + (73 cos a 3 t + 74 sin a 3 t)x + cos c^te 

1 • 
+7-^5 7r( a 2 cosa 3 t + a 2 sma 3 t 

1 



9 



3 3 sin a^t — /?4 cos a 3 t + —^5 ^- (a 2 sin a 3 t — a 3 cos c^^x — 

(CKj + QJg) / at 



+ — <5 4 cos a 3 t + <5 3 sin a 3 t + (73 sin a 3 t — 74 cos a 3 t)x + sin a 3 tx 2 



X 7 



1 • ^ A d 

+ , — — («2 smc^t — a 3 cos a 3 t)x — — 



[a\ + a 2 3 ) 



dx 



1 \ d 

fa cos a 3 t + /3 4 sin a 3 t + 7— 2 ^( a 3 sm a ^ ~ a ? cos ot^t)x ) 77: 



dt 



+ ^5 cos a 3 t + ^6 sin a 3 t + (75 cos a^t + 76 sin a 3 £)a; + cos a 3 tx 2 



X* 



1 • ^ A d 

+ ; — — - (a 3 sm a 3 t - a 2 cos a 3 t)i J — 



(a 2 . + a 2 ) 



1 \ d 

/?4 cos a 3 t — /?5 sin a 3 t + 7— ^ -37 ( a 2 sin a 3 t + a 3 cos a 3 t)x ) — 



(«2 + «3 



+ 5q cos a 3 t — £5 sin a 3 t + (7 6 cos a 3 t — 75 sin a 3 t)x — sin a 3 ta 2 



1 • x A 9 

+ / o — tt ( a 2 sm a 3 t + a 3 cos a 3 t)x — — 



(«2 + «1) 



<9x 



(59) 



As far the authors knowledge goes the above symmetry generators are being reported for 
the first time. Eq. ( 1551) is also linearizable one as it admits eight parameter Lie point 
symmetries. The general solution^ of Eq. (|55|) reads as 



, . —3a 6 
x(t) = + 



(3a 2 + A 1 )(l-/ 1 e ±ia *) 



k k V 3(1 - J ie ±«) - 2(3a 2 + Ai)/ 2 e(- 3a± f >* ± d(l - he ±&t ) 



(60) 



where 



3 + a\\ — k\ 2 = 0, and a = \/— 3a 2 — 4Ai. 



(61) 



For more details about the derivation of the solution (160]) one may refer to Ref. 
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B. Generalized Modified Emden Equations: / = k\x + k 2 



Next we consider the case where / as the most general linear function of x, that is 
/ = k\X + k 2 , where k\ and k 2 are arbitrary constants. As a consequence we arrive at 

9 = + + \ x x + A 2 , (62) 

where Ai and A 2 are arbitrary constants. 

Proceeding in a similar fashion as earlier, we consider the following four cases: (i) X\ = 
X 2 = 0, (ii) Ai 7^ 0, A2 = (Hi) X 2 7^ 0, Ai = 0, and (iv) Ai 7^ A 2 7^ 0, separately. The 
corresponding equations of motion, respectively, are 

x + (k l x + k 2 )x + ^^x 2 + ^-x 3 = 0. (63) 
3 9 

u k 9 k 2 

x + (k x x + k 2 )x + -^x 2 + -±x 3 + X x x = 0. (64) 
3 9 

x + (hx + k 2 )x + ^x 2 + fx 3 + A 2 = 0. (65) 
3 9 

fci ko k 2 

x + (k x x + k 2 )x H x 2 + —x 3 + Aix + A 2 = 0. (66) 

3 9 

However, all the above equations can be transformed straightforwardly to the equations 
discussed in the previous subsections. For example with the subtitution, X — x + jj^, Eq. 
( 163]) can be transformed to the form 

P P 9 P 

X + hXX + ^A 3 - ^-X + ^ = 0. (67) 
9 3 9 k\ 

Equation (1671) coincides exactly with (|55j) by restricting the constants Ai = — f and X 2 = 
|p in the latter. As a result the infinitesimal symmetries and their associated symmetry 
generators can also be derived from the Eqs. (1581) and ( |59l) . 
Eq. (164j) can be cast into the form, 

i HlX x4^ + !^X + l^W=0, (68) 

through the transformation A = x + Eq. (1681) is of the form (1551) . Therefore, the 
infinitesimal symmetries and the associated symmetry generators of equation (|68|) can also 
be derived from the Eqs. ( |58i) and ( |59i) by appropriately fixing the constants. 
In a similar way, one can transform Eqs. ( 1651) and ( 1661) to the following form: 



X + hXX + § .Y 3 - M.Y + M + ^ = 0. (69) 
* + + ^ *3 + (A, - |)X + ^ + 9k j - 9hM = 0. (76) 
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respectively. Eqs. (|69|) and (1701) are also of the form (1551) and hence we can derive the 
infinitesimal symmetries directly from the Eq. fl58l) . 



V. NON-MAXIMAL LIE POINT SYMMETRIES FOR f xx £ 
1. General forms of f(x) 

One can consider higher degree polynomials or other forms for / such that (f xx ^ 0) 
and classify the equations and their underlying Lie point symmetries. However, our analysis 
shows that one can get only lesser number of Lie point symmetries than eight whenever 
fxx 7^ 0. This can be proved in a rather general way as given below. 

Now we assume / = k±x + k 2 + l(x), where l(x) is an arbitrary function and investigate 
the outcome when Z(x) ^ 0. With this choice, g takes the form (vide Eq. (jSJ)) 
1 



9= 3bl 



ix 1 - 2bf 2 + cx + d)f x - (a + bx) + 46/i - 3bx - 2c + a . (71) 



In the above we introduced the notation fx — J f(x)dx and f 2 — f f\(x)dx— J (J f{x)dx)dx. 
Substituting g and its derivative into Eq. (j4j) and expanding the latter we arrive at an 
expression exclusively in terms of l(x), that is 

4,/, k 2 n k x n\ 2 c 2 lx 2 2cdlx Id 2 4 / A fk 2 9 ki » 
bl(U + -fx 2 + -rx 3 — + - Z cx + d)[ -fx 2 + -rx 3 + Z 4 



3 \ 2 6 J 3b 3b 3b 3 V / V 2 6 

+ ~3~b + + — 3~ (^ 2X + + ^ 3 ) + + ^ 2 ) _ ~b + ^) (° X + / 
+2 (y x 2 + y x 3 + Z 4 ) {k x + /i) a - ^ (& 2 + fcix + Z 2 ) - (fox + y x 2 + Z 3 ^l - 



+ jj& (fox + y x 2 + Z 3 ) 2 + ^ (fo + fox + Z 2 ) (4cx + 3d) - (fo + Zi) - ^ (fo + Zi 

-4(yX 2 + yX 3 + Z 4 ) (fo + fox + Z 2 )fe - jjx(fox + yX 2 + Z 3 ) (fo + fox + Z 2 ) 6 ~ ^ 

+2a; f ^x 2 + ^x 3 + Z 4 ) f A* + 6 - -bx 2 l (cx + d) + -x 2 lb { ^x 2 + ^x 3 + h) - ^ 
V2 6 /V / 3 \ / 3 V 2 6 J b 

+ yr(fox + yx 2 + Z 3 j - y|-(fo + kxx + Z 2 )x - y (fo + Zijx 2 + y (fo + fox + Z 2 jx 2 

6 2 x 3 /, ,\ Z6 2 x 4 2cc 4./, fo 2 , \ . /, 

; — fo + Zi : 1 : c /c 2 X H x + Z 3 + cx fc 2 + fox + Z 2 

6 \ V 36 36 3 V 2 / V 

/ \ ; cd 2 / fo 2 \ 2dd cZfe •■ / 2 \ 

+ (^fo + fox + Z 2 Ja - — + -(^fox + yx +^ 3 J a + ~^ ^ — 26(^fox + fox + Z 3 Jx 

+x 2 (k 2 + fcix + Z 2 ) fc i + cx + d + bx 2 = 0. (72) 
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Again in (1721) . for the sake of simplicity, we have used the notation l\{x) = fl(x)dx, 
— J [J l(x)dx]dx and so on. It is clear from the above expression that the first term 
on the left hand side of Eq. (1721) . namely |6/^/ 4 + ^-x 2 + ^-x 3 ^j cannot be balanced by no 
other term. This in turn leads us to 

^/(/ 4 + yX 2 + |x 3 ) 2 ^6 = 0. (73) 

When l(x) ^ the above expression clearly indicates that the symmetry function b(t) should 
be zero. Since one of the symmetry functions becomes zero one can get ultimately lesser Lie 
point symmetries only. This has in fact been analyzed in detail in Sec. Ill of paper I. 

On the other hand restricting l(x) = 0, that is / = k\x + then one obtains maximal 
Lie point symmetries as shown in the earlier sections. 

2. Relationship between the symmetries and the forms of f(x) and g(x) 

In our discussions so far, we find that all the examples dealt with in Sees. Ill and IV 
admit eight symmetry generators, and each one of the sets form sl(3,R) algebra. It is a 
well known fact that second order ODEs with eight symmetry generators are linearizable by 
point transformation^^, which was pointed out by Lie himself^. In fact the theorem that 
provides criteria for testing the linearization of a scalar second order ODE 'x = f(t, x, x) via 
a point transformation has the cubic in derivative forrn^^ 

x + A(t, x)x 3 + B(t, x)x 2 + C(t, x)x + D(t, x) = 0, (74) 

with the coefficients A,B,C and D satisfying the following two invariant conditions, 

3A a + 3A t C - 3A X D + 3AC t + C xx - 6AD X + BC X - 2BB t - 2B tx = 0, (75) 
6A t D - 3B X D + 3AD t + B tt - 2C tx - 3BD X + 3D XX + 2CC X - CB t = 0, (76) 

where suffixes refer to partial derivatives. On comparing Eq. (0Q) with Eq. (1741) . we find that 
the functions A(t,x) = B(t,x) = and C(t, x) and D(t,x) are f(x) and g(x) respectively 
which are independent of time. Substituting these forms of C(t,x) and D(t,x) along with 
A(t,x) = B(t,x) = in Eqs. £3} and (J76D, we obtain 

fxx = 0, 3g xx - 2ff x = 0. (77) 
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Integrating Eq. (1771) . we obtain 

/ = k\x + k 2 , g = -hk 2 x 3 + -k\x 2 + \\X + A 2 , (78) 
y o 

where Ai,A2 are integration constants. Thus for the function f xx ^ 0, one obtains lesser 
parameter symmetries only . However, our analysis leads to a more general form from which 
one can determine systems exhibiting lesser parameter symmetry groups. Thus this analysis 
broadens our knowledge on both integrable and linearizable equations belonging to Eq. (CD). 

Note that the above relation (1781) is the same as the one obtained from the analysis of 
the determining equations, Eq. (21) for the Lie symmetries in Sec. II A of paper I for the 
case fxx ~~ 0* 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied in some detail the symmetry properties of the Lienard 
type equation (1), where none of the symmetry functions is equal to zero. We have also 
isolated the class of equations which possess eight Lie point symmetries. For the sake of 
completeness we have also given the solutions for all the nonlinear ODEs which we have 
identified as linearizable in this procedure. We have also indicated the explicit Hamiltonian 
forms for these equations. 

A question now naturally arises is about the Noether symmetries of these equations. As 
it is well known, Noether's theorem is based upon the invariance of action integral and 
essentially one needs a Lagrangian to start the analysis. Since our recent studies show that 
even dissipative systems can have time independent Lagrangian 21 and so it is worth exploring 
Noether symmetries of these systems as well. Finally, it is also of importance to study the 
detailed symmetry structure when the functions / and g in (1) are also time dependent. 
These are being pursued currently. 
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APPENDIX A: METHOD OF SOLVING THE DETERMINING EQS. (f50j) 



In this section we solve the determining equations (joOl) to obtain the general form of £ 
and i] which leaves Eq. ( 1491) invariant. The determining equations to be solved are 

36 + ka + kc = 0, (Al) 

2c - a + 3A 2 6 + kd = 0, (A2) 

c + kd=0, (A3) 

d - A 2 c + 2A 2 a = 0. (A4) 

In the following, we solve eqs. flAlfl - flA4J) with the assumption that the parameters k, A 2 7^ 
and deduce the infinitesimal symmetries. 
Rewriting (I All) we get 

3 •■ 

c=--b-d. (A5) 
k 

Substituting (IA5I) into (IA2I) and rearranging, we obtain an expression for the function d of 
the form 

t 6 ••• 3 .. 3A 2 / . _\ 

d= ¥ b + k a ~~ k (A6) 

Now substituting the forms of c and d and their derivatives into (1A3I) and (1A4j) and simpli- 
fying, we get two equations for the unknowns a and b: 

d A b , <i 3 a , ^ db , „ s 

3 dF + 2k ^- 3kx ^t =0 ' (A7) 

<i 5 6 , <i 4 a , , da , . s 

6 rf? + 3 *dF + 3 ^S = a < A8) 



From Eq. flA7j) we have 



<i 3 a 3 d 4 6 3 db . . 

^3- = ~2l^ + 2 A2 ^- (A9) 
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Substituting flA9|) into (lA8|) we get 



d 5 b o , _ d 2 & da 

Rearranging the above we obtain 

da 1 d 5 6 3 d 2 b 

dt ~ ~2k 2 X 2 d^ ~ 2k d^ ' 

On integration (lAlip yields 



a(t) = a\ 



1 d% 3 d6 



2fc 2 A 2 dt A 2k dt ' 

where a± is an integration constant. Using (lAllj) . Eq. (1A7j) can be written as 
A general solution for (IA13I) takes the form 



b(t) = b 1 + + + ^ + ^ + ^e»* + ^e»*, 

Pi P2 P3 Pi P5 Pa 



(A10) 



(All) 



(A12) 



(A13) 



(A14) 



where b^s, i = 1, . . . , 7, are integration constants and p'jS, j = 1, . . . , 6, are given by 

Pl = (-3)^(a 2 )5, p2 = -(-3)s(A:A 2 )5, p 3 = (-l)i(-3)^(a 2 )3, (A15) 
p 4 = -(-l)*(-3)s(A;A 2 )s, p 5 = (-l) § (-3)s(A;A 2 )3, p 6 = -(-l)§(_3)s(fcA 2 )s. 

Substituting the form of 6 given by f)A14|) . and its derivatives into (IA12I) . we get 

I 3); 



a(t) =a x - {3+ u ] b2e^ - { - 
(3 + v^3) 



2k 



2k 



6 4 e P3t 



hePi t _ (3 + ^3) ^ _ (3 V3) 67eP6 t_ 
2 2 /c 2 k 

Using (1A14|) and (1A16I) . an explicit form for c can be derived from the relation (1A5j) : 



(A16) 



c(t) 



-3 + 



-3) 



+ 



2A; 
-3 + 



pib 2 e 



pit 



(3 



-3) 



2A; 



p 2 b 3 e 



Pit 



(3 + v^) 



2k 
(3 4 



=3) 



p 6 6 7 e P6 *. (A17) 



2 /u 2 2/c 

Now with the known expressions for the functions a, b and c one can deduce the form of d(t) 
through the relation ( 1A6I) . namely, 



d(t) 



™2 o, n «t , m l 



2k 2P5 6 2k\ 6 ' 



(A18) 
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where 



(3 + v^3) (3-^=3) /Ai . 

mx = 2fe2 ; , m 2 = ^ } . (A19) 



For further analysis we redefine the constants 



7712 9 m l 2 m l 2 

ZjIxi-^ Zifv-^ Lti\j^ 



m 2 2 m 2 2 m l 2 

r 4 = ^pP4, r 5 = —2 P5 , r 6 = — 2 P6 . (A20) 

£rv-^ £ihj^ Zirv^ 



and write the function d(t) in the form 

d(t) = r^&i + nb 2 e Plt + r 2 b 3 e P2t + r 3 he Pst + r 4 b 5 e P4t + r 5 b 6 e P5t + r 6 6 7 e P6 *. (A21) 

k 



APPENDIX B: METHOD OF SOLVING THE DETERMINING EQ. flggh 

Here we solve the equations 

36 + kd + 3A X 6 + kc = 0, (Bl) 

2c-d + 3A 2 o + A;d = 0, (B2) 

c + 2Aid + kd = 0, (B3) 

d-A 2 c + 2A 2 d + A 1 d = (B4) 

k 2 

b +k\ 2 b + X 1 b + kc+ —d = 0, (B5) 

and deduce the general form of a(t), b(t), c(t) and d(t). 
From ([Bl]) and we get 

3- 3 

c = ——b — —Xib — a, (B6) 
k k 

6 ••• 6 , 3 .. 3A2 , /T ^„ N 

d = v b+ ¥ Xlb+ k a -— b - m 

By substituting (1B6[) and (1B7I) and their derivatives into (1B3[) and (1B4h we obtain 

d 4 6 , d 2 & , , db , d 3 a , , da , 

3 5? + 3A ^- 3M2 5 + 2 ^ + 2M '*= ' (B8) 

d 5 6 d 3 fe „,,(!() , d 4 a , , d 2 a . , 2 da rt 

2 ^ + 4Al ^ + 2A? * + + **' + w * = a (B9) 



From (1B8[) we have 



d 3 a da 3 d 4 5 3Ai d 2 b 3\ 2 db 

d£3 + 1 dt = ~2^dF ~ ~2k~dl? + ~~dl' ( ^ 
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Substitute flBTUj) into flB9]) and rearranging we get 

da 1 d 5 b 5Ai d 3 b 3 d 2 b 2A 2 db 

dt ~ ~2k 2 \ 2 dtf ~ 2k 2 \ 2 dtf ~ 2k d^ ~ l?\ 2 ~dt 

Integration of (IB 111) leads to 



(Bll) 



1 d% 5Ai d 2 b 3 db _ 2Xj 
a[t) ~ ai ~ W\ 2 d£ ~ W\ 2 W ~ 2k dt ~ ¥\ 2 ' { ' 



where a\ is an integration constant. Substituting (1B1 1|) and its derivatives into (1B8|) we 
arrive at 

d 7 b , d 5 b ^, 2 d 3 b _, 2 . 2 ,(I6 n 

^ + 6Ai 5F + 9a; 5F + (4A ' + * a3 * =0 - (Bl3) 

Introducing 4| = P, where P is the new dependent variable, we can reduce an order of the 
Eq. (jBT3l) . that is, 

^ + 6A ^ + 9A ^ + < 4A ? + 3fc2A ^ < Bl4 > 

Eq. (IB14I) is a sixth order linear ODE with constant coefficients whose solution can be found 
in the following way. Since Eq. ( ]Bl4j) contains constant coefficients a general solution of this 



equation can be written of the form 

P = p ie mit + P 2 e m2t + P 3 e m3t + P 4 e mit + P 5 e mst + P 6 e m6 *, (Bl5) 

where Pi, i = l,...,6, are integration constants and m^s, i = l,...,6, are roots of the 
characteristic equation 

m 6 + 6Aim 4 + 9A 2 m 2 + (4A? + 3k 2 \ 2 2 ) = 0. (Bl6) 



Now choosing m 2 = X we can rewrite the sixth power polynomial equation (IB16I) as a cubic 
polynomial equation in X, namely, 

X 3 + 6A1X 2 + 9A 2 X + (4A 3 + 3k 2 X 2 2 ) = 0. (Bl7) 

Let us choose 

a = 6Ai, 6 = 9A 2 , c = A\l + 3k 2 \ 2 2 , (Bl8) 
and introduce a transformation 22 

X = Y ~l ( Bl9 ) 
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in (IB 171) so that it becomes 



Y 3 + pY + q = 0, 



(B20) 



where the new constants, p and q are related to the old constants, a, b and c, by the following 
relations 

,2 



a' 2 o ab 

p = b , q = — a he. 

y 3 ' ^ 27 3 



(B21) 



The cubic equation (IB20I) has one real root and two complex conjugate roots^, namely, 

Yl = A + B, K„ = Jd±2 ± ,^Z*>, (B22) 

where 



•4= (-| + a/Q) S - « 



9 



2 ■ «=&i nt 



p 



(B23) 



so that the roots of Eq. (IB 1 Tj) can now be expressed through the relation (jB19p . that is, 

(A + B) a .y/3(A-B) 



X 1 = A + B--, X 2 



- + i- 
2 3 



(A + B) a .VS(A-B) 

v\q — — — — I 

2 3 2 



(B24) 



From the identity m 2 = X, we find the roots of the characteristic equation (IB16I) : 



a 



m lj2 = ±y / X 1 = ± X /A + B-~, m 3A = ±^X 2 = ±VC + iD, 



±VX 3 = ±VC - iD, 



(B25) 



where again for simplicity we have introduced the constants C and D which can be fixed 
from the relation 



c _ (A + B) a _ V3(A - B) 
2 3' 2 



(B26) 



Eq. (1B25j) can be put in the following convenient way 



mi 

m 3 
m 5 



■JA + B-^, m 2 = -^A + B-^, 

(C 2 + D 2 )^ [cos0 + isin0] , m 4 = -(C 2 + D 2 )^ [cos0 + 2sin</>] , 
(C 2 + D 2 ) L * [cos0 - zsin0] , m 6 = -(C 2 + D 2 )^ [cos0 - isin</>] , 
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where 



tan" 



(B27) 



2 J 

With the help of these m^s, i = 1, . . . , 6, a general solution for (1B14I) can be written in the 
form 

P = b 2 e ait + b 3 e~ ait + e a2t (fo 4 cos a 3 t + b 5 sin a 3 t) + e~ a2t (b 6 cos a 3 t - b 7 sin a 3 t) , (B28) 

where b' { s, i = 2, . . . , 7, are integration constants and a 2 , a 3 are real and imaginary parts of 
m 3 ,m4,m 5 and m 6 respectively, where a 2 = VC 2 + D 2 cos (ft, 03 = VC" 2 + D 2 smcj). 
From the identity ^ = P, the function upon integration can be deduced as 

b 2 I) 646°^^ 

bit) = bi H e ai * e~ Ql * + —^5 srfaa cos a 3 t + a 3 sin a 3 t) 

ai q;i (a 2 + 0:3) 

6 5 e a2 * b e e~ a2t 
2 2T ( a2 s ^ n a ^ ~ 0:3 cos a3 ^) + 2T ( a 3 sm a ^ 



3> 



— a 2 cos a 3 t) + =- (a 2 sin a 3 t + a 3 cos a 3 t) , (B29) 

(«2 + "3) 

where 61 is an integration constant. Once b(t) is known the function a(t) can be fixed from 
the relation flBT2l . 



2A 2 

a{t) = ^ - — + fab 2 e ait + fab 3 e~ ait + b±e a2t {fa cos a 3 t + fasma 3 t) 

K X 2 

+b 5 e a2t (fa sin a 3 t — fa cos a 3 t) + b 6 e~ a2t (fa cos a 3 t + /3 4 sin a 3 t) 

+b 7 e~ a2t {fa cos a 3 t - fa sin a 3 t) , (B30) 

where a\ is an integration constant. For simplicity in the expression a(t) we have introduced 
the constants fa, i = 1,...,5. These new constants are related to the old parameters 
through the following relations 



, oj 5Aiai 3_ _2A|\ = / a\ hMai _ 3_ 2A 2 

~" ~ \2PA 2 2k 2 X 2 2k k 2 X 2 aJ ' ^ 2 V 2 fc 2 A 2 2fc 2 A 2 2/c k 2 X 2 ai 

/a 2 (a2-3a§) 5Ai« 2 3 2A 2 a 2 

— — ^TTT r 7T7^ r- 



2fc 2 A 2 2FA 2 2k k 2 X 2 (a 2 2 + a 2 3 

,'a 3 (al - 3a 2 ,) 5Ai« 3 2A 2 a 3 

— — 77T9\ r- 



2k 2 X 2 2k 2 X 2 k 2 X 2 {a 2 + a 2 ) 

_ f a 2 (a 2 2 - 3a 2 3 ) 5Ai« 2 3 2A 2 « 2 , { 

Pb ~ V 2FA 2 2fc 2 A 2 2A; + PA 2 (a 2 + a 2 ) ' ' 
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With the forms of a(t) and b(i) the function c(t) can be deduced using the relation (1B6I) in 
the form, 

3A 

c(t) = + 7i& 2 e ai * + l2he~ ait + 6 4 e a2 *(7 3 cos a 3 t + 74 sin a 3 t) 

K 

+65 e° 2 * (73 sin a 3 i — 74 cos a 3 t) + b§e~ a2t (75 cos a^i + 76 sin a 3 £) 
+6 7 e _a2< (7 6 cosa 3 /f: — 7 5 sino; 3 t), (B32) 

where again we have defined new constants ^s, i = 1, . . . , 5, for simplicity sake. Here too 
one can fix their forms using the following relations, 



7i 



3«2 3AiQ;2 _ _ \ / 3a3 3AiQ!3 



73 = - ~r + i 7 2 I 2T + "2/^3 + a 3 /3 4 , 74 = -J ? / 2 , 2T + "3$? ~ "2^4 

3a2 3Aia2 „ „ \ / 3a 3 3AiQ!3 



75 = + ^fT^I) +a2/?5 ~ a3/?4 J' 76= l~~ fc(^ + ^) +a2/?4 + a3/? 7' 

Finally, the function d(t) can be derived from the relation flBT|) by simply substituting the 
forms of a(t) and b(t) and their derivatives into it, that is, 
3A 

d(t) = r^foi + 5i6 2 e ai * + 5 2 b 3 e- ait + he a2t (8 3 cos a 3 t + 5 4 sin a 3 t) 

+b^e a2t {5 3 sin a 3 t — £4 cos a 3 t) + 6 6 e~ a2 * (5s cos a 3 t + 56 sin a 3 t) 
+b 7 e~ a2t (5 6 cos a 3 t - 5 5 sina 3 t). (B34) 

Again for simplicity we have defined the constants S^s, i = 1, . . . , 6, which can be determined 
from the following relations, 



1 / Qa\ 6A1 2 3A 2 \ . If 6af 6A1 2/3 3A- 



h = T ( 3 («2 - «3)(^5 + 7) + ^ - 6«2«3A + 7 3 2 Q2A2 2 

re \ k k [a.\ + a 



3/ 



5e = 7 ( 6a 2 a3(^5 + 7) + 3/3 4 (a2 - of) - J • ( B35 ) 

re \ re ^2 + a 3 J / 
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